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1 Introduction 



As is well-known, the Adler-Bell-Jackiw anomaly Jll in the divergence of the axial 
vector current in QED, and in the axial flavor singlet current in QCD, is intimately 
connected with the fact that a shift of the integration variable in the linearly divergent 
Feynman integral for the triangle graph does not leave the integral invariant. An 
additional surface term is generated by the shift which, if one requires vector current 
conservation, leads to an anomalous axial vector Ward identity. By introducing a 
Pauli-Villars regularization such shifts in momentum variables are allowed, but chiral 
symmetry is explicitly broken by the Pauli-Villars regulator masses. Alternatively one 
can use the dimensional regularization scheme of 't Hooft and Veltman to compute 
the anomalous contribution to the divergence of the axial vector current. 

An alternative regularization scheme is that provided by the lattice formulation of 
QED and QCD. Using Wilson fermions one avoids the fermion doubling problem 
at the expense of breaking chiral symmetry explicitly, as demanded by the Nielson- 
Ninomiya theorem |Q. The term in the action which is responsible for lifting the 
fermion degeneracy is a so called irrelevant term which vanishes in the naive continuum 
limit. Since the lattice regularization is already introduced on the level of the generating 
functional, the lattice Ward identities for the divergence of the vector and axial vector 
currents are expected to hold for any finite lattice spacing. While it follows that the 
divergence of the vector current is conserved, the divergence of the axial vector current 
is not. The corresponding lattice Ward identity yields an explicit expression for the 
anomalous contribution, which is entirely determined from the Wilson term in the 
action that breaks the chiral symmetry explicitly. This anomalous contribution, which 
is given by an irrelevant operator, has been calculated a long time ago by Karsten and 
Smit 1^ and was shown to yield the correct anomaly in the continuum limit. The axial 
anomaly on the lattice has been discussed subsequently by various authors 0. 

In the continuum formulation of QED or QCD the anomaly must necessarily be 
computed from the divergence of the axial vector current, since the naive Ward identity 
does not provide us with an explicit expression, as is the case on the lattice. One is then 
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led to the conclusion that computing the dimensionally regulated triangle graph for 
the divergence of the axial vector current in the continuum formulation is equivalent to 
taking the continuum limit of a lattice Feynman integral having no continuum analog. 
The fact that the anomaly computed from an irrelevant lattice Feynman integral yields 
the triangle anomaly of the continuum formulation is not obvious. There is no a priori 
reason why this should be the case. It is the purpose of this paper to show why this is 
necessarily so. 

The paper is organized as follows. In section 2 we consider the axial vector Ward 
identity for the U (1) lattice gauge theory with Wilson fermions, and check explicitly 
that the anomalous contribution to the divergence of the axial vector current is given 
by an irrelevant lattice integral which has precisely the form demanded by the Ward 
identity. An important ingredient in the proof is that lattice Feynman integrals are 
invariant under shifts in integration variables. In the continuum formulation, on the 
other hand, such shifts in linearly divergent Feynman integrals generate an additional 
surface term which give rise to the anomaly, as we have already remarked before. It is 
this observation which provides us with the key for proving that the continuum limit 
of the anomalous contribution to the Ward identity is necessarily given by the triangle 
graphs in the continuum formulation. The basic idea is that, having proved that the 
Ward identity is satisfied for any finite lattice spacing, we can also write the irrelevant 
anomalous term in the lattice Ward identity as the sum of contributions of the indi- 
vidual Feynman diagrams to the divergence of the axial vector current without making 
a shift of integration variables. This sum will now include, besides the additional dia- 
grams arising from irrelevant lattice vertices, the lattice version of the familiar triangle 
diagrams which do possess a naive continuum limit. In section 3 we then apply the 
small-a expansion scheme of Ref. to this expression and show that for vanishing 
lattice spacing it reduces to the dimensionally regulated integrals for the triangle di- 
agrams of the continuum formulation. In order to point out a similarity between the 
mechanism by which the anomaly is generated on the lattice and in the continuum 
formulation, we shall apply the small-a expansion scheme also to the alternative ex- 
pression given by the irrelevant anomalous term in the lattice Ward identity. Section 
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4 summarizes our results. 



2 Check of the lattice Ward identity 

Consider the action for lattice QED with Wilson fermions 

4r^ 



S 



Wilson 
F 



1 



[ipj^j] = ( m H j'^tp (x) tfj [X) 



-— J2 (^) - 7m) Uf, (x) i/j {x + afi) + ip {x + ajj,) (r + 7^) [/^ (x) ip (x) 



, (2.1) 



where J2x = J2n ?n is the fermion mass, r the Wilson parameter and a the lattice 
spacing. From the invariance of the partition function under the change of variables 
ip (x) — > exp[iA (x) 75] tp (x), and ip (x) ^ ip (x) exp[iA (x) 75] one readily derives the 
following Ward identity for the expectation value of the divergence of the axial vector 
current in an external field {U^ (x)}, 



< A^jj^^ (x) >=2m <ip (x) -f^ip (x) > + < X (x) >, 



(2.2a) 



where is the left lattice derivative. The axial vector current is given by 

-^i^^ (a;) = 2 ^ (^) 7M75f^M {x) ^p {x + afi) + ^p {x + afi) ^i^^i^Ul (x) ^ {x) 
and the irrelevant operator X (x) by 



(2.2b) 



X (x) = — — ^ 4^ (a^) 75 U^, (x) -0 (x + afi) + [/^ (x — a/i) tp {x — aji) — 2^p (x) 
Za ^ 



— — ^ (x — a/t) f/^ (x — a/i) + ip {x + ajl) Uj^ (x) — 2ip (x) 
za ^ L 

On the other hand, the vector Ward identity reads 



75^ (x) . (2.2c) 



<Ajj^(x) >=0, 



(2.3a) 



where 



Jf^ (x) = 2 



^ (x + afi) (r + 7^) Lr[ (x) ^ (x) - ^ (x) (r - 7^) (x) (x + a/i) 



. (2.3b) 



Since the lattice regularization is introduced (non-perturbatively) at the level of the 
partition function, the Ward identities (|2.2a| ) and ( |2.3a|) should hold for any finite 



lattice spacing. If so, one can calculate the anomalous contribution to the divergence 
of the axial vector current directly from < X (x) >, as was done originally in Ref. P|. 
This contribution is proportional to the Wilson parameter r, with X (x) an irrelevant 
operator that vanishes in the naive continuum limit. The explicit computation of 
< X (x) > turns out to yield the correct anomaly for vanishing lattice spacing P . 

In the continuum formulation, on the other hand, the anomalous Ward identity 
must necessarily be computed from < d^jjj^^ (x) > with jj^^ (x) = ifj (x) (x). We 

shall follow here the same procedure, and verify in one-loop order that this external 
field expectation value is in fact given by the RHS of eq. (|2.2a|) for any finite lattice 
spacing and Wilson parameter r. The purpose of this exercise is to show that the proof, 
that the anomalous part of < jj^^ (x) > is actually given by an irrelevant Feynman 
integral, rests heavily on the fact that lattice regulated Feynman integrals are invariant 
under shifts in integration variables (which is not the case, e.g., for linearly divergent 
continuum Feynman integrals). It is this observation which will provide us with the 
key to establish the equality between the anomaly calculated from the dimensionally 
regulated triangle graph in the continuum formulation, and the anomaly obtained from 
the expectation value of the irrelevant operator X (x) for vanishing lattice spacing. At 
the same time our check of the lattice Ward identity gives a nice example of the role 
played by irrelevant vertices in ensuring eq. ( |2.2a| ). 

Consider first the LHS of the Ward identity (g^. By setting [/^ (x) = e^^"^"(^+'^A/2) 
in eq. ( |2.2b| ), and expanding this expression up to second order in the potentials 
(notice that they are defined at the midpoints of the links) one is led in one-loop order 
to the Feynman diagrams depicted in Fig. 1. The fermion propagator on the lattice 
reads 

(q) =I^Miq)+'-Y^ 7, sin q,a^ , (2.4a) 
with M (q) the momentum dependent Wilson mass, 

M{q) =m + -Y (1- cos qpa) . (2.4b) 

a y 

The lattice expressions for the vertices coupling one or two photons to the vector 
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current are given by 



Vx VP) = (7a cos — ir sm — ) , (2.5a) 



(p) = le'aS,, (7A sin ^ + cos ^) , (2.5b) 

where p = p + p', with p [p') the incoming (outgoing) fermion momentum at the 
vertex. In the hterature it is customary to write the external field expectation value of 
A^J^^-* (x) in the form 

- — ^ sin ^——^Ti^xu {k, k'; a) = fxu {k, k'; a) , (2.6) 

where k and k' are the outgoing photon momenta, and T^xu is computed with the 
following expressions for the axial vector vertices, denoted by dots in Fig. 1, involving 
no, one or two photons: 

4m (P) = 7m75 cos ^ , (2. 7a) 

4m,a (p) = -«e(5^A7M75 sin (2.7b) 



(p) = -^a^e2(5^A5;,^7M75 cos (2.7c) 

For the following discussion it turns out however to be important to consider the 
corresponding axial vector vertices for Txu written in a particular convenient way. 
These can be obtained directly from the definition of the current (|2.2b|) : 



{p,p; a) =ia 7^75 (sinp^a - sinp^a) . (2.8a) 



Wx^ (P, P, k; a) = ie-fx^o cos {2p -k)^^- cos {2p' + k)^ ^ . (2.8b) 



^A? (P'P'i k, k'; a) = -^ae%A75^Ai. sin {2p - k - k')^ | - sin i2p' + k + k')^ |j (2.8c) 
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It can be easily verified, making use of trigonometric relations, that the lattice Feynman 
integrals corresponding to diagrams (a)-(f) in Fig. 1, computed with these expressions, 
can be written in the form ( p.6|) with the effective axial vertices (|2.7a|) - (|2.7c| ) for F^aj/- 
Consider now the contribution of diagrams (f) and (c). Clearly diagram (f) does 
not contribute to the divergence of J^^ (x). Diagram (c) can also be easily shown to 
vanish by making use of the following identity for (|2.8a|) 



la 



'E7M75(sinp^a-sinp;a) = -75D(p)-D(p')75 + 75[M(p) + M(p')]. (2.9) 



Next consider the contributions arising from diagrams (a), (b), (d) and (e). Making 
again use of the identity (|2.9|) we can write the contribution of diagrams (a) and (b) 
in the form 

Ti^+') (A;, k'; a) = [ -^N,, {i, k, k'; a) + Mx, {k, k'; a) , (2.10a) 



with 

Nx. {e,k,k';a) 



^^(l5Vi^^ (2£ - k) D-^ {£) Vj^^ (2£ + k') {£ + k') 
+Tr (^75^-1 {i - k) Vi^^ {2i - k) D-^ (£) KJ^^ {2^ + k') ) 
+ [{k,\)^{k\v)l (2.10b) 



and 



M\y {k, k'; a) = - 



M{i-k) + M{i + k') 



L (2-)' 

a 

xTr(^75D~^ (£ - A;) y^f^^ (2£ - k) {^) V^^^ {21 + k') D'^ + k!) ^ 
+ [(A;,A) ^ (A;',z/)]. (2.10c) 

Here [(fc, A) ^ (A;', v)\ denotes the contribution of diagram (b). Since the term propor- 
tional to the Wilson parameter in the Vertex v'^'^ (p) [cf. eq ( |2.5a| )] does not contribute 
to the trace in N\y (£, k, k'\ a), we can make the following replacement in eq. ( |2.10b| ) 



Vl'\p)^-ie[jx + nxm], (2.11a) 
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where 



vanishes in the continuum hmit. Hence we obtain 



(2.11b) 



Nx, {e,k,k';a) = 

= i-ief [ti (^75[7a + m - k) ]D-' {£) [7, + {2i + k') ]D-' {£ + k') ) 

+Tr (75^-^ (£ - k) [7A + ^^A (2^ - k) ]D-' (£) [7, + fi, (2£ + k') ]) 
+ [{k,\)^{k\v)]. (2.12) 

This expression can be decomposed into irrelevant terms (i.e., vanishing in the contin- 
uum hmit) proportional to one or two Vt^ (9)'s, and terms which possess a continuum 
analog. By making appropriate shifts of integration variables these latter terms can be 
readily shown not to contribute to the integral ( |2.10a| ). In the continuum formulation 
the very same shifts would also lead naively to a vanishing result, but because one 
is faced with linearly divergent integrals a surface term is generated which is nothing 
else but the anomaly. In a lattice regulated theory such shifts of variables leave the 
integrals invariant and the anomaly is actually generated only by M.\y {k, k'; a), given 
in eq. ( p.lOc|) . In fact, as we now show, the contributions of the irrelevant terms in 
eq. ( p.l2| ) to the integral ( |2.10a| ) are canceled by the contributions arising from the 
lattice diagrams (d) and (e) involving the irrelevant vertex w^^^^ {p,p',k;a) defined in 
eq. (|2.8b|) . Writing w^^^ {p,p', k; a) in the form 

{P,p',k;a) = -ie75 



Qx {2p -k)-Qx {2p' + k) 



(2.13) 



where VLx {q) has been defined in eq. ( |2.11b| ), and noting that we can make again the 
replacement ( |2.11aD for the vertex coupling the photon to the vector current (the term 
proportional to r in eq. (|2.5a]) does not contribute to the trace), we obtain the following 
expression for the contributions of diagrams (d) and (e) 



Xv 



{k, k'; «) = / Tr|^75[r]^ {2i - k') - {2i + 2k + k') ]D-^ (£) 

a 

X [7A + nx (2£ + k) ]D-^ (£ + A;) )1 + [(A;, A) ^ (A;', v)]. (2.14) 



By appropriate shifts of integration variables one then readily verifies that the contri- 
bution of these diagrams cancel the remaining irrelevant contributions in eq. (p.l2|) 



to the integral in ( p.lOa| ). Hence the LHS of the Ward identity (|2.2a| ) is given just by 



Aixu, defined in eq. (|2.10c| ), i.e. by an irrelevant lattice integral. It is now easy to 
check that this result for < Aj^jj^^ (x) > coincides with that computed from the RHS 
of the Ward identity ([Oal) . Since the vertices associated with the irrelevant operator 



X (x), defined in eq. ( p.2c|) , are proportional to 75, only diagrams (a) and (b) will now 
contribute because of the symmetry properties of the Dirac trace. The corresponding 
75-vertex is given by 



Mr{p)+Mr{p') 



(2.15a) 



v{p,p;a) = 75 
where 

7W^(g) = -^(1 -cos^po) (2.15b) 
a p 

is the momentum dependent Wilson mass ( ^.4bD evaluated for m = 0. Including the 
contribution of 2m < -0 {x) -j^ifj (x) >, the Feynman integral corresponding to the RHS 
of eq. ( p.2a| ) is therefore given by (|2.10c|) . 



Notice that only after using the identity ( |2.9| ), and making appropriate shifts in 
integration variables, we were led to the conclusion that the contribution of diagrams 
(a)-(f) to the anomalous part of < A^jj^^ (x) > combined to yield the expectation 
value of the irrelevant operator X (x) . On the other hand, the divergence of the axial 
vector current written in the form (|2.6| ) with the effective axial vertices ( p.7aD -( p.7c|) 



for T^xu, does not only include the irrelevant contributions of the diagrams (c)-(f), but 
also the contributions of diagrams (a) and (b) which do possess a continuum analog. 
Since we have shown the Ward identity to be satisfied for any finite lattice spacing, we 
have therefore the following relation 

fxAk,k';a)=MxAk,k';a), (2.16) 

where Txu {k, k'; a) and M. xu (k , k' ; a) have been defined in eqs. (|2.6|) and (|2.10c| ). 



respectively. In the following section we will now use the small-a expansion scheme of 
Ref. [0 to show that in the continuum limit Txu reduces to the dimensionally regulated 
Feynman integral for the triangle diagrams of the continuum formulation. 
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3 Computation of the axial anomaly in the small-a 
expansion scheme 

In the small-a expansion scheme of Ref. one departs from a D-dimensionally regu- 
lated one-loop lattice Feynman integral having the generic form 

+- 

F {p,m;a,D) = / ^ H {i,p,m; a, D) , (3.1) 

a 

where p and m stand collectively for the set of external momenta and masses, and 
it is understood that the coupling constant e is to be replaced, as is usual in the 
dimensional regularization scheme, by e/x*-^"^^/^, with /i an arbitrary mass scale. An 
important feature of the lattice Feynman integral is that the lattice spacing appears in 
the integrand as well as in the integration limits. This integrand is a periodic function 
of p and p' and reduces to the integrand of the corresponding continuum Feynman 
integral in the limit of vanishing lattice spacing. 

According to the small-a expansion scheme F [p, m; a, D) has the following power 
series expansion in the lattice spacing:^ 



+ 00 

F {p, m; a, D = 4) = lim 



D^4 



f d^i 

-oo 

+a'^--^T,_(,^_4) J H (i, ap, am; a = 1, d) ] + O (a^+(^-^)+i) , (3.2) 



where du is the inverse-mass dimension of H {i,p, m; a, D), i is the dimensionless loop 
variable, and Tkf (a) denotes the Taylor expansion of / (a) around a = up to O (a*')- 
Note that the integration in the first integral extends, as in the continuum, over an infi- 
nite range, and that the leading O (a°) contribution is just the dimensionally regulated 
continuum Feynman integral. The Taylor expansion of the second integral yields poly- 
nomials in the momenta with coefficients having the form of lattice Feynman integrals 
with periodic integrands. 

In the following we will first use the small-a expansion scheme to calculate the 
anomaly from eq. ( |2.10c| ). Although this computation is not required for the main 

*The small-a expansion has been applied also in [Q, 
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objective of the paper, it serves to illustrate the distinct ways in which the relevant 
operator 2m'?/' (x) ^^ij) (x), and the irrelevant opeiaXoi X (x) in the Ward identity (|]2D 
are handled by the small-a expansion scheme. At the same time it will also serve to 
point out the similarity between the mechanism by which the anomaly is generated on 
the lattice and in the continuum formulation. 

i) The anomaly calculated from the small-a expansion of Aixi, {k, k'; a) 



Consider the RHS of the Ward identity ( p.2aD . As we have seen, it is given in one- 
loop order by Ai\u {k, k'; a), defined in eq. ( ^.lOcI ). Extracting from this expression the 
contribution proportional to m, making use of eqs. (|2.4b|) and ( p.l5b|) , we have that 



Ai\u {k, k'; a) = 2m 



where 



7D/ 



(2vr) 



D 



I^£{i,k,k';a) + I^^^{e,k,k';a) 



Axu (k, k'; a) , (3.3a) 



Axu {k, k'; a) ■- 
+- 

d^i 



(2n 



.D 



Mr {i-k)+ Mr + k') ] I^^J {i, k, k'- a) 



+ [{k,X)^ik\u)i{3.3b) 



is the anomalous contribution to the divergence of the axial vector current. Here 

{i, k, k'; a) = 

= -Tr (^-f^D-^ {£ - k) Vi^^ {2£ - k) D'^ (£) Vj^^ (2£ + k') D'^ {£ + A;') ) , 

and 

lSii,k,k';a) = ll''^{i,k',k;a). (3.3c) 

The fermion propagator D"^ (q) and the fermion-photon vertex V^^^ (p) are defined 
in eqs. (|2.4a|) and ( p.5a| ), respectively. We now apply the small-a expansion scheme 
( |3.2| ) to this expression, with dn = 2 in the present case. Since it turns out that the 
corresponding Taylor coefficients in ( p.2|) , up to O (a°), are finite for D — 4, we can 
set D = 4 in this equation. 

Consider the integral in ( p.3a|) whose leading contribution to the first, as well as to 
the second term in (R.2|) is found to be of O(a^). The contribution in this order to 
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the second term vanishes due to the symmetry properties of the Dirac trace. We are 
therefore only left in O (a°) with the contribution of the first term in eq. ( p.2|) , which 
just yields the usual expression for the divergence of the axial vector current arising 
from the explicit breaking of chiral symmetry due to the fermion mass m|]] 

= 2mTxv {k, k') , (3.4a) 

a 

where 

2 1 ^ 



lim 2m 



^xJ (^5 k,k 



(£, k, k ; 



What concerns the anomalous term Axi, [k, k , a), in eq. (3.3a), exactly the opposite 
is true. Since the integrand in ( p.3b|) vanishes in the continuum limit, it does not 
contribute in O (a°) to the first integral in eq. ( p.2| ). It does however contribute to the 
second term in this equation. Although the leading term in the expansion is of O (a^^) 
one finds, after a straightforward, but lengthy calculation, that the coefficients of 
and vanish, and that 

+7r 



A\,y {k, k') = lim ^ (fc, k'; a) = 16e^ exu,ipkf,k'p 



(3.5a) 



where 



Hx, (I) 



cos ix COS 4 COS if^iXif. (^tj COS ip - ArM,. {Pj sin^ tp 



(3.5b) 



and where all quantities with "hat" are measured in lattice units. This coincides with 
the expression of Karsten and Smit 0. After making use of the identity 



(3.6) 



Ml {i) COS ip - ArMr (i) sin^ ip = cos ip [M^ [i) +Asm'^ i 

+ (mi {i) + J2 sin' L) ' sin ^4- (m^ {i) + ^ sin^ 4 

^ K ^ Uin \ K 



Although the authors in Ref. j|] claim, without giving a proof, that they can also generate this 
contribution, we do not see how this is possible in their approach. In particular we question the 
remark following eq. (5.42) in their paper. 
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and performing a partial integration one readily verifies that 

i 

V-Vt 

where 



A.(A:,A;')=lini / ^d,ICpxJlk,k') (3.7a) 



^p\u A;') = lQe^exy^,pk^k'p cos 4 cos 4 cos 4t 



smCp 



2 ' 



(3.7b) 



where the integration extends over a 4-dimensional cube with edges lying in the first 
Brillouine zone [ — tt, +7r], excluding an infinitesimal sphere with radius e around the 
infrared singularity at the origin.^ By applying Gauss theorem and taking the limit 
e — > 0, one obtains the result of Ref. p 

(fc, k') = -^sx^ppk^k'p, (3.8) 

which is just the well-known anomaly of the continuum formulation. This, as we have 
already mentioned before, is a remarkable result, since the anomaly has been computed 
from an expression having no continuum analog. On the other hand, the sum of the two 
triangle diagrams contributing to the anomaly in continuum QED are known to give a 
non-vanishing contribution only because linearly divergent Feynman integrals are not 
invariant under shifts in integration variables. The surface term generated by such 
shifts, arising from an ultraviolet divergence, yields the anomaly (see e.g. Ref. 0). 
On the lattice the surface term is generated by a similar mechanism, but the anomaly 
arises from an infrared divergence. This, and the above explicit computation, suggests 
that one can actually prove that the anomaly arising from the explicit breaking of chiral 
symmetry by the irrelevant Wilson term in the action is necessarily identical with that 
calculated from the continuum triangle Feynman diagrams. The crucial ingredient 
going into the proof is, that since we have shown that the Ward identity ( p.2a| ) is 
satisfied for any finite lattice spacing, we can just as well compute the anomaly from 
the divergence of the axial vector current, using the small-a expansion scheme without 



Since the integral in ) is infrared finite, we are fi-ee to calculate it as the limit e of the 



volume integral V—Ve- This provides a natural regularization of the, in the limit e ^ 0, logarithmically 
divergent integrals appearing after making use of eq. ( |3.6|) . 
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making use of the identity ( p.9| ) and of the invariance of lattice Feynman integrals 
under shifts in integration variables, i.e. from eq. ( p.l6|) , where Tx^ has been defined 



in eq. (|2.6| ). Since the Taylor expansion of the individual lattice Feynman integrals 
contributing to Txu are sensitive to shifts in integration variables, we will be able to 
transfer the anomaly arising from an irrelevant lattice operator, to the familiar triangle 
anomaly in the continuum. We now show how the small-a expansion scheme is able to 
realize this program. 

ii) Alternative computation of the anomaly 

Consider the individual Feynman integrals associated with the diagrams (a)-(f) con- 
tributing to Txiy in eq. ( p.6|) with the momenta chosen in the way indicated in the Fig. 
1. Since the integrands of the integrals associated with diagrams (a) and (b) possess 
a continuum limit, while those associated with the remaining diagrams vanish in the 
limit a — » (due to the irrelevant vertices), we immediately conclude that the leading 
contribution of O (a°) of the first term in ( ^.2| ) just yields the usual dimensionally reg- 
ularized continuum expression for the triangle graphs (a) and (b). The computation of 
the contributions to the second term in this equation, which involves only lattice inte- 
grals with periodic integrands, is very tedious. The coefficients of O (a^^) and O (a^^) 
turn out to be finite for D A, and vanish because of the symmetry properties of the 
trace. In O (a°), however, one is also faced with infrared divergent lattice integrals for 
D — s> 4 at intermediate stages of the calculation. Hence one must consider the dimen- 
sionally regulated expressions, as demanded by the small-a expansion scheme ( p.2|) . In 
the 't Hooft-Veltman scheme p[ {"ffj.,'ju} = ^df^u remains valid in D dimensions, while 
75 anticommutes with 7^ for fi = 1, ■ ■ ■ , 4, and commutes with 7^ for > 4. In the case 
where one is dealing with convergent lattice integrals for D —>■ 4, one can work directly 
in four dimensions using the usual anticommutation rules for the 7^- and 75-matrices. 
What concerns the in the limit D ^ 4 divergent integrals, we found, after a lengthy 
computation, involving many partial integrations, and making only use of the anticom- 
mutation relations {7^, j^} = 26^^, that one is left with the following expression for 
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the Taylor coefficient for the contribution to the second term in eq. ( |3.2| ) in O 



(3.9a) 



where 



V, 



1 COS COS Ia cos^ 


-3Ml 




p 


(2nf 


Ml (i) 


+ E sin^ ip 
p 


3 



(3.9b) 



and 



Q 



COS COS ix COS 


4 sin^ ix - E sin^ ^ 
p 


i2.f 


^ ^ p 


3 



(3.9c) 



Since the traces in eq. (p.9a|) vanish for D — 4, the only non-vanishing contribution 
to this coefficient can arise from the divergent parts of the integrals (|3.9b|) and (|3.9c| ). 
Consider ffist the expression ( p.9b| ). The potentially divergent contribution to this 
integral arises from the integration region {£p} ^ 0, where cosf^ ~ 1. What concerns 
sin^^K, it is convenient to write it in the form sin^i^ ^ 4sin^ y. With these approx- 
imations the integral in eq. ( p.9b| ) becomes independent of A, so that we can further 
replace sin^ ^ by E sin^ y. Then the potentially divergent contribution reduces to 



D 



4 7 d'^l 



D J (2n 

— TT ^ 



4Esin2| 



2 • 



For D — i> 4 the divergent part of this integral behaves like [Stt^ {D — 4)]^^ (cf. Ref. |^). 
We therefore find that the integral (|3.9b|) is actually convergent for D — 4, so that 
the ffist term in ( p.9a| ) can be evaluated directly in four dimensions and vanishes since 
Tr (757^/7^) = for D = 4. Proceeding in a similar way one finds that the contribution 
of the second term in eq. (|3.9a| ) also vanishes for D —>■ 4. Hence, in this alternative 



computation of the anomaly, the second term in eq. (3^) does not contribute in the 
continuum limit, and the anomaly, as well as the contribution (3.4) arising from explicit 
breaking of the chiral symmetry by the fermion mass, are obtained entirely from the 
ffist integral in eq. (p.2|). 
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We have therefore shown that taking the continuum hmit of the anomalous contri- 
bution to the divergence of the axial vector current, arising from the irrelevant operator 
X {x) in the Ward identity (|2.2a| ), is necessarily equivalent to taking the D —>■ A limit 
of the dimensionally regulated continuum Feynman integrals for the triangle diagrams. 
Hence the small-a expansion scheme has allowed us to establish a direct connection 
between the anomaly calculated from an irrelevant lattice Feynman integral and a 
dimensionally regulated continuum integral. 

4 Conclusion 

In this paper we have taken a new look at the axial anomaly in lattice QED with 
Wilson fermions, originally discussed in Ref. 0. We have first checked explicitly that 
the axial Ward identity is satisfied in one-loop order for any finite lattice spacing and 
Wilson parameter. In checking the lattice Ward identity, we have made extensive use 
of the fact that lattice Feynman integrals are invariant under a shift of integration 
variables, which is not the case in general for ultraviolet divergent integrals in the 
continuum formulation. By making use of this property we have shown explicitly that 
the anomalous contribution to the divergence of the axial vector current is given in 
one-loop order by an irrelevant lattice Feynman integral. At the same time this check 
provided a nice example of the role played by the irrelevant vertices in ensuring the 
validity of the Ward identity. Having shown that it holds for any finite lattice spacing 
we were then led to consider an alternative expression for the anomalous part of the 
divergence of the axial vector current. By applying the small-a expansion scheme of 
Ref. to this expression we have then shown in a systematic way that the axial 
anomaly, as computed from an irrelevant lattice operator, necessarily coincides with 
that computed from the dimensionally regulated triangle diagram in the continuum 
formulation. The crucial ingredient going into the proof was that the small-a expansion 
of lattice Feynman integrals is sensitive to shifts in integration variables. This made 
it possible to transfer the anomaly arising from an irrelevant lattice operator to the 
triangle anomaly computed in the continuum formulation of QED in the dimensional 
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regularization scheme. 
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Figure caption 

Figure 1: Lattice diagrams contributing to the divergence of the axial vector current. 



17 




(d) (e) (f) 

Figure 1: Lattice diagrams contributing to the divergence of the axial vector current. 
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